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Abstract 

We show that Brans-Dicke (BD) theory in 5D may explain the present cosmic accelerated expansion without 
recurring to matter fields in 5D or dark energy in 4_D. Without making any assumption on the nature of the 
extra coordinate or the matter content in 5D, here we demonstrate that the vacuum BD field equations in 5D are 
equivalent, on every hypersurface orthogonal to the extra dimension, to a BD theory in 4D with a self interacting 
potential and an effective matter field. The potential is not introduced by hand, instead the reduction procedure 
provides an expression that determines its shape up to a constant of integration. It also establishes the explicit 
formulae for the effective matter in 4D. In the context of FRW cosmologies, we show that the reduced BD theory 
gives rise to models for accelerated expansion of a matter-dominated universe which are consistent with current 
observations and with a decelerating radiation-dominated epoch. 
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1 Introduction 



One of the most challenging problems in cosmology today is to explain the observed late-time accelerated expansion of 
the universe [J- [13]. Since the gravity of both baryonic (ordinary) matter and radiation is attractive, an accelerated 
expansion requires the presence of a new form of matter (called dark energy), which could (i) produce gravitational 
repulsion, i.e., violate the strong energy condition; (ii) account for 70% of the total content of the universe; (iii) 
remain unclustered on all scales where gravitational clustering of ordinary matter is seen. (For a recent review see 
Ref. [II]). Possible candidates for dark energy include: a cosmological constant or a time dependent cosmological 
term [15j-[16j: an evolving scalar field known as quintessence (Q-matter) with a potential giving rise to negative 
pressure at the present epoch [T7]-[1D]; dissipative fluids [^T]; Chaplygin gas [11]-[13; K-essence [11]-[17], and other 
more exotic models [55]. 

General Relativity (GR) is a well tested theory on solar system scales. In contrast it is poorly tested on cosmic 
scales [29^. Therefore, the question arises of whether dark energy is not just an observational artifact caused by 
an inappropriate theory of gravity. This and other puzzles of theoretical and experimental gravity have triggered 
a huge interest in alternative theories of gravity (See, e.g. [20] and references therein). In recent years there has 
been a renewed interest in scalar-tensor theories of gravity as viable alternatives to general relativity. In particular, 
some researchers [31]- [37] have resorted to Brans-Dicke theory (BD) [38] in order to explain the present accelerated 
expansion of the universe. An attractive feature of BD is that the scalar field is a fundamental element of the theory, 
as opposed to other models in which the scalar field is postulated separately in an ad hoc fashion. The concept 
is that the BD scalar field could play the role of Q-matter or iiT-essence and lead to cosmological acceleration. 
However, it turns out that this is so only in very particular cases: for values of the coupling parameter w in the 
range —2 < uj < —3/2, which not only violate the energy condition on the scalar field but are also inconsistent with 
a radiation-dominated epoch, unless uj varies with time j32| ; when there is a scalar potential, which in turn is added 
by hand [33]. Further models include the so-called "chameleon fields" that allow the scalar field to interact with 
matter (see [39] and references therein). 

Another attempt to derive the present accelerated expansion from a fundamental concept (or postulate) combines 
the original BD and Kaluza-Klein (KK) theories with the modern view that our AD universe can be recovered on a 
hypersurface orthogonal to the extra dimension. Along these lines, and assuming that the extra coordinate is cyclic 
(ignorable), it has been shown in [40] that the BD field equations in 5D are equivalent to those of GR in AD with 
two scalar fields, viz., the BD scalar field and the 755 component of the KK metric (For a recent discussion see Ref. 
[41) ). These two scalar fields may account for the present accelerated expansion of the universe, if one admits the 
existence of matter in 51? which does not move along the extra dimension and effectively behaves as dust in AD. 

In this work we adhere to the point of view advanced in [40]- [41]. However, we propose that neither the cylindricity 
condition on the extra coordinate nor the higher dimensional matter hypothesis are necessary. In fact, in this paper 
we show that, regardless of whether we assume a compact or large extra dimension, BD theory in 5D may explain the 
present accelerated expansion without introducing matter fields in 5D or dark energy in AD. We demonstrate that 
the vacuum BD field equations in 51? are equivalent, on every hypersurface orthogonal to the extra dimension, to a 
BD theory in AD with a self interacting potential and an effective matter field. In the context of FRW cosmologies, 
we show that the reduced BD theory gives rise to models for accelerated expansion of a matter-dominated universe 
which are consistent with a decelerating radiation-dominated epoch, for the same value of the BD parameter a;. 

The paper is organized as follows. In section 2 we perform the dimensional reduction of the scalar-vacuum BD 
field equations in 5D. Our procedure provides explicit definitions for the effective matter and potential in AD. In 
section 3 we study homogeneous and isotropic solutions of the vacuum BD field equations in 5D under the assumption 
of separation of variables (in general, the metric functions may depend on time and the extra coordinate). We show 
that there are four classes of separable solutions. We discuss in detail the class of power-law solutions, which in 
turn consists of five families of solutions. In section 4 we study the AD cosmological scenarios allowed by the 5D 
power-law solutions derived in section 3. We find that distinct families of solutions in 5D lead to very much alike 
scenarios in AD; they can give accelerating matter-dominated era as well as decelerating radiation-dominated era for 
the same value of uj in the range — 3/2<ti;<— 1. In section 5 we give a summary of our results. Finally, in the 

'^A similar procedure has recently been discussed in 1421 . where a "modified" BD theory in 4D is obtained from the vacuum BD 
equations in 513. However, our formulae and definitions of the AD quantities are completely different from those derived in | 42| . 



2 



Appendix we show, by means of explicit integration, that the assumption of separabihty is consistent with the field 
equations. 



2 Dimensional reduction of Brans-Dicke theory in 5D 

The Brans-Dicke (BD) theory of gravity in 5D is described by the action [3D] 



(5) = / d'x^/\^\ 



1671 / d''x^^\Lg\ 



(1) 



where i?^^-* is the curvature scalar associated with the hD metric 7ab; t''^' is the determinant of ^ab] </> is a scalar 
field; w is a dimensionless coupling constant; Lm'' represents the Lagrangian of the matter fields in bD and does not 
depend on (f). 

The equations for the gravitational field in 51? derived from ([1]) read 



^AB - ^AB " -jlABrL^ - —J A R 



AB + ZH 



(Va0) (Vb0) - --fAB (V^0) (Vc</ 



(VaVb0-7asV2,^) , (2) 



where = V^V"* and Tj^g represents the energy-momentum tensor (EMT) of matter fields in 5-D with trace 
y(5) — j^^T^g. For the sake of generality of the reduction procedure in this section T^^ ^ 0. 
The field equation for the scalar field is determined by ([T]) as 



2a;. 



Taking the trace of ^ we find 



Combining the last two equations we get 



A + 3oj 

In this work we use coordinates where the metric in 51? can be written a^ 

dS^ = ^Asdx^dx^ — gp,u{x, y)dx^dx'^ + e$^(x, y)dy'^, 



(3) 

(4) 
(5) 

(6) 



in such a way that our AD spacetime can be recovered by going onto a hypersurface "Ey : y — yo — constant, which 
is orthogonal to the 51? unit vector 



-A 04 



— , UAU =e. 



(7) 



along the extra dimension, and 17^1, can be interpreted as the metric of the spacetime. 

The effective field equations (FE) in 42? are obtained from dimensional reduction of ([2]) and ([5]). To achieve such 
a reduction we note that 



e * 
2*2 



e 

$2 



2 $ 



(8) 



■^Notation: a;'' = (x'', x^, x^, x^) are the coordinates in 4D and y is the coordinate along the extra dimension. We use spacetime 
signature — , — , — ), while e = ±1 allows for spacelike or timelike extra dimension. 
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where the asterisk denotes partial derivative with respect to the extra coordinate (i.e., d/dy = *); Da is the covariant 
derivative on Ey, which is calculated with g^i,, and = D°'Da. 

Using these expressions, the spacetime components {A — = v) of the bD field equations ^ can be written 

as 



Stt 



T(5) + — 



2$2 



(9) 



To construct the Einstein tensor in AD we have to express and i?*^^^ in terms of the corresponding AD 

quantities. The Ricci tensor r[^} of the metric 5^,/ and the scalar field $ are related to the Ricci tensor i?^^ of jab 
by i43j 



- 9 al3 +9 9aX9p^ 



44 4 2 

From ©-(HI) and the second equation in (ITUl) we obtain 



A;3 



2$ 



(10) 



2$2 



$0 



A/j ** , 9 9x0 $ g^'^ 5a/3 
9 5 A/3 H 



$20 



* * 

** * I LO (j) $ 



4 + 3w $2 
Substituting this expression into R'^^^ — j'^^Rab we find 



(11) 



2 (i^q$) e 

4$2 



5 3a/3 + (5"^ 9afij 



2e 







/ * * \ - 








- 


f 





167r 



(12) 



$2 4 + 3a; 

where i?'-^^ = g^^R^^p is the scalar curvature of the spacetime hypersurfaces S^. 

• We are now ready to obtain the effective equations for gravity in AD. With this aim we substitute the first 
equation in (ITUl) and into © and isolate g|5 = R^f} — g^^R^^^^ /2. The result can be written as 



(D^^) [D,<i>) - ^5^, {Dac^) (D- 



i {D^D,^ - g^.D^<f>) - 5m. (13) 



where we have introduced the quantity V^('/>), which (as we will see bellow) plays the role of an effective or induced 
scalar potential; S^i, is the reduced EMT of the matter fields in 5D 



4 + 3w 



44 
$2 



(14) 
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and T^a^^^ can be interpreted as an induced EMT for an effective BD theory in AD. It is given by 



with 



(15) 



^rpiSTM) ^ 



2$ 



9 9ap + [g"'^ 9af}j 



(16) 



Since in BD acts as the inverse of the Newtonian gravitational constant G, (fT6|) is identical to the induced EMT 
used in STM (Space-Time-Matter theory) [43 . The second term in (fT5|) depends on the first derivatives of with 
respect to the fifth coordinate and represents the effective EMT in AD coming from the scalar field. The equation 
for the scalar potential V is given bellow by ([19]). 

Taking the trace of ^3]) we obtain a simple relation between i?^'*), S = g'^^S'^^ and T^^^) = 9'"'tI^^\ namely 
(we note that g'^'Tin) = T^^' - e^ilV*^) 



V 



3D^ 



2V 



(17) 



To construct the AD counterpart of ([5]) we substitute (fT^ and ([T7| into ([3]). After some manipulations we get 



3 + 2w 



1 



3 + 2w 



21/(0) 



where 



(18) 



2(l+w) 



^ $2 



e0 

4$2 



5 .9a/3 + (ff"'^ 5a/3 




(19) 



This equation, with the r.h.s. evaluated at some 'Sy, constitutes a working definition for the potential. 

We notice that equations (|13p and ([T5)) are identical to those of Brans-Dicke theory in AD derived from the action 



(4) 



IGtt / d x^^gL 



(4) 
m 7 



(20) 



with 5 



5, 



(BD) 



2dSm/Sg'^'^ , where Sm — J d'^x^J—gh'^hi represents the action for matter in AD. 



The potential vanishes only in few particular cases. For example when cj = — 1, y is a cyclic coordinate and 
= 0. However, if we are to recover a general version of BD theory in AD, we must assume uj ^ —\. Besides, the 
metric in bD can depend on y. Therefore, in general V ^ Q. 

• Finally, for the line element ^ we have Rfl = $Pf.. [13] with 



1 

2$ 



(21) 



The dynamical equation for P^p is obtained from ^ by setting A — A, B = viz.. 



D,. 



(22) 
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In the case where T^4 = and — constant, this reduces to P".^ — 0. In braneworld theory this quantity is 
proportional to the EMT of the matter on the brane If ?/ is a cychc coordinate, (1^^ reduces to = 0. 

To keep contact with other works in the hterature, we note that we can substitute (ITT|) in (ITil) to obtain an 
expression where the right hand side of (1131) only contains T^li/ as well as the derivatives of 0, $ and g^jy. In case that 
there is no y dependence we neatly recover the formulae developed in [40 . However, our equations are completely 
different from those derived/used in [42] , 



3 Scalar- vacuum Brans-Dicke cosmology in 

In cosmological applications the 51? metric ([6]) is commonly taken in the form 



1 — kr^ 



{dO^ 



e^\t,y)dy^ 



(23) 



where k = 0,+l,— 1 and {t,r,9,(j)) are the usual coordinates for a spacetime with spherically symmetric spatial 
sectionfU. 

For this line element the vacuum (74^ = 0) Brans-Dicke field equations ^ reduce as follows: The temporal 
component A = B = gives 



,d / a $ 

5- - + 



3fcn2 ien^ 
a? $2 



a a I a $ 
a a \ a $ 



1 




$2 



n 2(f) 



(24) 



the spatial components A = B = 1,2,3 reduce to 



2d 


d ( a 


2n\ 




$ /2a 


fi 


a 


h- -- 
a \ a 


n J 




'A a- 


n 



the A ^ B ~ 4 component gives 

3fcn2 



d 

- H 




(d 


--)\ 


a 


a 


\a 





the mixed component A = 0, B = 4 yields 



kr? 



Finally, the wave equation ([S]) becomes 



en 
"$2" 



3en I a \ I a n 
$2 \ a I \ a n 



2 a a I a 2 n 
a a I a n 



n $ / 2 a n 
n $ I a n 






n (f) 
ruj) 



LjJ 

la' 



en (j> I a oj 











* / a; ^ 


-:)] 









(26) 



(27) 



V^0 = — 




e 

$2 




(28) 



^We do not make any assumption about e. Rather we let the FE to determine the signature of the extra dimension. In Cases 1 and 5 
(with the exception of solution II68I 1) the extra dimension must be spacelike. All other solutions allow both signatures. 
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3.1 Separation of variables 

In this section, we look for sohitions to the above equations under the assumption that the metric coefficients are 
separable functions of their arguments. In this framework, without loss of generality we can se10 

n{t,y)=N(y), a{t,y) ^ P{y)Q{t), ^t,y) ^ F{t), ^{t,y) ^ U{y)W{t). (29) 

Then, from (|27|) we obtain 

0. (30) 



3Q 


(k 


-i) 


•fi 




-i) 


w ( 


Q ' 


Ijv 






{ p ' 







u 



Thus, separability requires 



N P_ f SP U_ f N {l+uj)U _ f 

N-p-^iJ^ — + ^ jj -C3-, (31) 

where ci, C2, cs are arbitrary constants and / = f{y) is some function to be determined by the field equations. 
From (|3T|) we get 



where A''o, Pq and Uq are constants of integration. Besides, from (pO| and ([3T|) it follows that 

Q^'^^i^'^W^^ = constant. (33) 



Now, to obtain an equation for / we substitute ([29|1 and (|32p into the wave equation (|28|) . The requirement of 
separability yields a second-order differential equation for /(y) whose first integral is 

^ (^Csf^"^ + C) j2[(3-3ci-4c2)w-3ci-5c2-C3+4]/(3w+4) ^g^-j 

where Cg is a separation constant and C is a constant of integration. This equation admits exact integration, in 
terms of elementary functions, in several cases. For example, setting C = we get 



, y(3-+4)/[(3c,+C2)(. + l)+C3] ^ (3^^ + {0J + 1)+ C3 ^ 0, 

giv/cT,^ (3ci + C2)(.. + l) + C3=0, OO. ^"^^^ 



Similar power-law and exponential solutions can be obtained for = or C2 = 0. 

Next, we substitute the metric functions with ^ and ([M]), into the FE ([Mll-dlH]). We find that these 

equations have the following structure 

H,{t) + voCH,{t) [f{y)]-''' +k{3Lo + 4) Hs{t) [f{y)f'' = 0, (36) 

where ijq = c| (5 -I- 4a;) + 2c2 [c^ + 3ci (1-1- a;)] — 3 (cs — ci)^; Hi{t) is a combination of F{t), Q{t) and W{t) and their 
derivatives; H2 oc W^Q^; cx W^F^. Inspection of the FE shows that the last two terms in (p6| do not cancel out 
for C2 = — ci 7^ and F oc Q. Therefore, the assumption of separability demands either ci = for /c ^ 0, or fc = 
for ci ^ as well as the fulfillment of (at least) one of the following conditions: C — 0, rjo — 0, C2 — 0. 

A detailed analysis of the field equations is provided in the Appendix, where we derive the differential equations 
to be satisfied by the three functions Q, F and W. Inspection of those equations reveals that they admit power-law 
and exponential solutions for several choices of the parameters. However, this is clear at the outset if we notice that 
(l30l) can also be written as 



*This is a consequence of the freedom to perform the coordinate transformation t = t{t), y = y{y). 
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which requires 



P [f Q j N [ Q W j U [f W 



F Q__h 3Q W__h F {l + uj)W 
F~Q~^^h' ~Q^W~'''^V F W 



= 0, 



(37) 
(38) 



where Ci, C2, C3 are arbitrary constants and h = h{t) is some function to be determined by the field equations. 
Integrating ([55)) we get 

p ^ f^l{3ci+C2){^+l)+C3]/{3uj+4) ^ Q /j[e2(w+l)-ei+g3]/(3c^+4)^ ^ ^[g2+3(ei-E3)]/(3tu+4) ^gg^ 

Following the same procedure as above, from the wave equation p8|) we find 



^ (C's/l^'^^ + C) /i2[(3-3Ei-4c2)w-3ci-5c2-C3+4]/(3w+4) 



(40) 



where and C are some new separation and integration constants, respectively. For several choices of the constants 
we obtain power-law and exponential solutions. For example, it we set C* = we get 



h{t) cx 



^(3^+4)/[(3ci+C2)(.. + l)+e3] ^ (3g^ +C2){OJ+1)+C3^0, 

(3ci + C2) (w + I) + C3 = 0, Cs > 0. 



(41) 



Once again, similar solutions are obtained = or C2 = 0. 

In summary, the assumption of separation of variables generates four distinct classes of solutions corresponding 
to the following choices: 



I 


hit) 


cx 


^(3w+4)/[(3ci+C2)(cJ+l) + 


^=•1, fiy) cx J/(3"+4)/[(3ci+C2)(l^+l)+C3 


II 


h{t) 


cx 


e±Vat^ /(2;)cxy(3 


^+4)/[(3ci+C2)(cJ + l) + C3] 


III 


hit) 


cx 


^(3w+4)/[(3ci+e2)(tJ+l)4 


^"^1, /(y)cxe±-^^ 


IV 


hit) 


cx 


e±^*, f{y) « e± 


5 



(42) 

where the constants obey the conditions indicated in ([55]) and (|1T|) . When we substitute (H^ in the FE the 
latter reduce to a system of algebraic equations which provide the appropriate consistency relations for the constants. 



3.2 Power-law solutions 

We now proceed to study in some detail the power-law solutions. We concentrate our attention on the spatially-flat 
scenario, which seem to be relevant to the present epoch of the universe j45j . 

We could use the analysis provided in the Appendix and study case by case the differential equations for F 
and W. However, since we already know the form of the desired solutions it is much easier to start from the family I 
in (j42p . For practical reasons it is convenient to adopt another parameterization, namely (01,02,03) — (a,/3,7) and 
(ci, C2, C3) — (/, m, s), where the power-law solution takes the form 

n^Ay"', a^ByH\ ^^Ct"", (jj^Dy^f. (43) 

Here A, B, C, D are some constants with the appropriate units; while a, /3, 7 and /, m, s are parameters that have to 
satisfy the field equations (l^- ipg)) . 
Substituting (|43| into (|28| wc obtain 

S (S - 1 + 3/ + m) C2t(^-2)y(T-2a) + _ ^ + 3^ + ^) ^2^(.-2m) ^(7-2) ^ q_ (44) 
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The above is satisfied in five different cases 



s 


= 0, 


7 = 


0, 




s 


= 0, 


7 = 


1 - 


3/3 -a, 


7 


= 0, 


s = 


1 - 


3/ — m, 


s 


= 1 - 


-31 - 


TO, 


7 = 1 - 3/3 - a, 


m 


= 1, 


a ~ 


1, 


s (s + 30 + £7 (7 + 3/3) ^2 



). (45) 

We now proceed to study these cases with some detail. 

3.2.1 Case 1: 

The selection s = 0, 7 = corresponds to the usual general relativity in 51? with — constant. In this case the 
power-law cosmological solutions are well-known in the literature j46j . However, to make the paper self-consistent 
we provide them here. 

From (|77|) we obtain al + /3 {m — I) ^ 0. Thus, assuming ^ 7^ to we get /3 = al/{l — to). In this case the only 
power-law solution to the field equations i^ 

dS^ = A^y^df - i32y2//(i-l)^2i ^^^2 ^ ^2 ^^g2 ^ gi^2 ^^^2)] _ -^^dy\ (46) 

(/ — 1) 

(The field equations cannot be satisfied for s = 0, 7 = and / = 1, simultaneously). For / — to the solution is 

dS^ = — dp - B^y [dr^ + (dO^ + sin^ ed(p^)] + edy^. (47) 

Here, contrary to what we see in P5|) . the coordinate y is not restricted to be spacelike. Therefore, by virtue of extra 
symmetry discussed in [47], the transformation of coordinates t — w^^^y, r = f/^/w, y = wt with w = ±A/e in (I47p 
generates the "new" solution 

dS"" = dp - BH [df^ + f2 {d9^ + sin2 ed^^)] + '-^df. (48) 

We should note that although the last two line elements are diffeomorphic in 5D (allowing complex transformations 
of coordinates in 51?), their interpretation in AD is quite differenl|f|. 

3.2.2 Case 2: 

The choice s — generates solutions where the scalar field depends only on y. Setting s — and 7 = (1 — 3/? — a) 
in p4)) we get m = 31 {a — f3) / {a — 1), where a 7^ 1. Now, p6| requires either I = 1/2 or Z = 0. 
The solution for I = 1/2 is 

dS^ = A^yf-'+'^y^dP - B^y^H [dr^ + {dO^ + sin^ Odip^)] + '-^dy\ 



^ = ,?,3a-5,)/4^ U+15.±,VW±M. (49) 

^ ' 94 + 75W ^ ' 

(When cj = — 94/75 « —1.25 the field equations require /3 = 3/5.) 

^We exclude here static universes where i = /3 = 0. These admit both signatures (e = ±1) and require one of the following: 
{0 = 1, m = 1); (0 = 0, m = 1); (a = 1, m = 0). 

®As discussed in section 2, in this paper we follow the approach where our AD spacetime is recovered by going onto a hypersurface 
'■ y = yo = constant, although other techniques are also possible 1481 . 
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The solution I — is 



dS^ = A'^y^°'dt^ - B'^y^'^ [dr^ + {dd^ + sin^ Bdip'^)\ + tdy^ , 

= , _lM!±i(5i^i)(f±i^. (50) 

(1 — a — 3/3) 

It should be noted that this solution works perfectly well for a — \ provided cj = —4/3. Also, when a = (1 — 3/3) the 
denominator in w vanishes and = constant. What this suggests is that the choice a — {\ ~ 3/3) should reproduce 5D 
general relativity. In fact, for this choice the field equations require either (a = —1/2, /3 — 1/2) or (a = 1, /3 = 0), 
which yield w — > oo. For the former set of values the Brans-Dicke solution (15(11) reduces to (H71) . while for the latter 
it reduces to the third static solution mentioned in footnote 5. 

3.2.3 Case 3: 

The choice 7 = generates solutions where the scalar field depends only on t. Setting 7 = and s = — 3/ — m) 
in (gll) we get a = -3/3 (l-m) j (m - 1), where m 7^ 1. From it follows that either ^ 1/2 or /3 = 0. The 
solution corresponding to /3 = 1 /2 is 

dS'' = ^de - B'^ye' [dr2 + r2 {dd'^ + sin^ 0^^] + eCH^^^^I^dy" , 

. Z^^3a-50/4^ ^^14+15.±4V6(4TM. (51) 
^ ' 94 + 75CJ ^ 

We note that this solution is related to (|49p by the transformation t ^ y, f3 ^ I, eC^ o A^. 
The solution for /3 = is 



dS^ = dt^ - S^i^' [dr^ + (rf6|2 + sin^ Odip^)] + eC^f'^dy^ 

(l-m- 30 



(52) 



which can formally be obtained from after the transformation t = (^^""C/A) y, r = w~^^f, y = wt with 
w = ±-y/e and replacing a and /3 by m and I, respectively. In the limit where m — )■ (1 — 3?), for which 4> = constant 
and cj — !■ CX3, we recover the general- relativistic solutions (I48p and the second static metric mentioned in footnote 5. 

It should be noted that ((52)) is the general power-law solution for the case where the metric depends on time 
alone. There is a very particular case, given bellow by (|57|) with /3 = 0, for which (j> is allowed to depend on y. 

3.2.4 Case 4 

We now look for solutions to the field equations of the form 



dS^ = A2y2(l-3;9-7)rft2 „ B^y2P^2l ^^^2 ^ ^2 (^^2 ^ ^-^^2 ^^^2)] ^ ^^2^2(l-3i-.) ^^2 ^ 

= Dy'^f. (53) 



Substituting these expressions into (|24p we get 

2 (72 - 7 + 37/3 -h 6/32 - 3/3) 



2 ,77^0 (54) 

7 



(If 7 = then we recover Case 3). Now from ((27l) we obtain 

_ 7^(1-60 + 37 [l + p{l-7l)] 
^~ 6/3 (1 - 2/3) 7 (7 + 1) ■ 



(55) 
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For these quantities we find that the choice /3 = —7/3 (s — —31) produces a simple particular solution, which is given 
bellow by (1551) . 

The general solution to the field equations in the case under consideration, is generated by a quadratic 

equation for I. The explicit expression for / is given by 



^_ 9(6/? + 77- l)/3^ + 3(8^ + 7- 2)7^-3(87 + 3) /3±V3r3/3 + 7) (3-6/3-7) [6 (2;3 - 1) ;g - (7 + 1) 7] , . 

3(473 - 1372 + 77-3) + 9[ll(6;3 + 57)/32-15(3;3-72)/3 + 4(2- 77)^] ' ^ ' 



Thus, (j54p -(|56 p describe a two-parameter family of solutions. From a mathematical viewpoint, this expression might 
look cumbersome and/or non-interesting. However, from a physical viewpoint, the fact that I explicitly depends on 
/? and 7 means that the expansion rate of the universe, which is determined by ^, is manifestly determined by the 
dynamics along the extra dimension. 

The above solution can be simplified if we add some physical requirements. The simplest ones arise from the 
choice of the coordinate/reference system. 

Synchronous reference system: The choice 500 = 1 is usual in cosmology; it corresponds to the so-called 
synchronous reference system where the coordinate t is the proper time at each point. Thus, setting A — 1 and 
a = (1 — 3/3 — 7) = in ([53]) the solution simplifies to 

dS^ = df - B^y^^t^' [dr^ + r2 {d9^ + sin^ 9d^^)] + eC^t^^^-^'-^'^ dy^ 

, . ny"-'^H: .^fi^. (57) 

where 



3 (3/3^ - 1) ± v/3 (2 ~ 3(3) [3/3 (/3 + 1) - 2] 

3 (3^2 + 6^-5) 
(1-3/3) [3/(l + /3)-l] 
3/3(l + /3)-2 



(58) 



It is worthwhile to mention that there are no solutions having simultaneously a = and /3 = 1/3. Besides, the 
condition /3 < 2/3 ensures that I is a real number and u > —4/3. For completeness, we note that there is a unique 
solution for which the metric depends only on time, while the scalar field is a function of both t and y. Such a 
solution is generated by the choice /3 = in the above expressions. 

Gaussian normal coordinate system: A popular choice in the literature is to use the five degrees of freedom 
to set (74^ — and 344 = 1. This is the so-called Gaussian normal coordinate system based on Ey. The solution in 
such coordinates can formally be obtained from ([57 1) - ([55)) by changing i o y, / O /3, s o 7; setting eC^ — A^; and 
allowing the extra coordinate to have one or the other signature. The result is 



61 {1-2 

(1 - 30 



^ ..^iii^, (59) 



where 



3 (3/2 - l) ± v/3 (2-30 [3^ (^ + 1) - 2] 



3 (3/2 + 6/ -5) 
3/(1 + /) 



7 . (60) 
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• Let us now go back to the general solution (1551) . For completeness we should consider the case where the 
denominator in (|55p vanishes. In such a case, from (l?7|) and ([55)) we find that the parameters /3 and 7 are restricted 
to be 

7 = -3, = 

^ 3/2 + 6^-1' 3^2 + 6/ -1' ^ ^ 

The first set of parameters generates two families of solutions. Namely, 

dS^ = A^y^df - B^y^f^ [dr^ + {d9^ + sin^ 9dip^)] + eC^i^^d,^) dy^, 

^ = Dy-h'i^^^\ "" = "1 (^2) 

where m(i) = 1, S(i) = —3/ and 771(2) = 3^ — 2, S(2) = 3 (1 — 21). For the second set of parameters in (j6ip . the expression 
for given by (15^ simplifies to a; = (1 — 9Z) /2Z and the field equations have solutions only for two particular values 
of /. These are I = 3/19 and I w 0.153, which is the real root of the cubic equation (246?^ - 171?^ + 40/ - 3) = 0. 

3.2.5 Case 5 

We now present the general solution for a = L m = 1. We assume s 7^ 0, 7 7^ 0, which were already considered in 
the previous cases. Thus, from (^5]) we obtairQ 

^2 £7(7 + 3/3) 2 , o, /Rox 



From (271) we get 

Now ([M)) requires either 



. = .. = ''<'-^)-^''="V'-^^"'-", (65) 

7S^ 



. = .2=^^^(^-^^ + ^^(^+;^^^^-^^r^^^^ + ^^(^ + ^\ 4. + 3-9/^0. (66) 

S7(4s + 3 + 90 ^ ^ 

Solution with lu — wi: Substituting (|55|) into ([M)) we obtain /3 = ^7/3. Consequently, from (p5)) it follows that 
= — e (7^/3^) yl^. Now, it can be verified that the field equations are satisfied by 

dS^ - A^y^dt' - B2y27//«i2/ ^^^2 ^ ^2 ^^^2 ^^^2 g^^2-)] _ 2!yl2^2^y2^ 

= Dy^t', Lo^uji. (67) 

The above solution works perfectly well for / = 1, although (IM)) requires Z ^ 1. Besides the extra coordinate has to 
be spacelike (e = —1). In addition, for the choice s — {I — 1)7 the line element (|67p becomes identical to the one in 
(H51) . but now (f) — Dy'^^^^~^h'^ and ui — (l — s)/s. With this choice, the 5D Brans-Dicke solution ([S7)) goes over the 
5D general- relativistic solution (|46p in the limil[f| a; — > 00 (s — > 0). 



^This case was recently discussed in 1491 . However, our solution 16711 seems to be much simpler. 

^To avoid misunderstandings, we should mention that the limit oj — >■ 00 does not necessarily implies that we recover general relativity. 
As an example consider the solution 1491 1: For ui —> 00 wc got /3 = 1/5 and (f> = constant, but the resulting metric does not satisfy the 
FE. A similar situation occurs in solutions 157I I- 158I I and 159I I- 160I I for /3 — )■ 1/3 and I — )■ 1/3, respectively. 
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For s — —31 the solution is 



4 



+ sin2 0V)] +eC'2t2dy2 



(68) 



We note that this is the only solution to the FE, with a = m = 1, which aUows the extra coordinate to be timelike 
(see (gni), dSZD). 



Solution with w — uj2: From the field equation ([^5]) we obtain the condition 

{s + 3lf[j{3l + l) + s{-f-l)] =0. 
Since s + 31 ^ (see dSS]) and ^) it follows that 

„ _ 7(3^ + 1) 



(69) 



7-1 



7 



1, / 



(70) 



It can be verified that when lu — uj2 the solution to the field equations is given by (|67p with the choice of parameters 
given above. When the denominator in vanishes, the solution is (|571) with 



3(1 + 3/) 



(71) 



Thus, the choice ui = uj2 does not generate new solutions to the FE. It just singles out some particular values for the 
parameters. Consequently, (p7)) is the general power-law solution with a = 1, m = 1. 

4 Reduced Brans-Dicke cosmology in 

We now proceed to study the effective AD picture generated by the 5D solutions discussed in the preceding section. 
For the line element (j23p . the cosmological metric induced on spacetime hypersurfaces T,y : y — y^ = constant is just 
the spacetime part of (|23p . and the no n- vanishing components of the induced Brans-Dicke energy-momentum tensor 
(fist are 



StT o(BZ5) 



Stt ,(BZ5) 



1 / $ n$ 













e 


n 








$2 


n 


a 




-1) 



2^' 



a a a * 2n 2 

- + :;: + — + 

a a \ a 9 n d 











n 


UJ(I)\ 






20 j 



V 

24,' 



(72) 



where V = V{(j)) should be determined from We note that r2(S£i) ^ ^3(515) ^ r^iiSD) ^ ^ffUch means that the 
induced EMT looks like a perfect fluid with energy density p = Tq'^^^^ and isotropic pressure p — 

Let us first calculate the potential. With this aim we substitute the power-law line element (|43| into (fT9|) and 



evaluate on Sj,. We obtain 



where 



-2ms [W + 1) i?2/^^-2j^2(7-as)/.^(,_2)/. _ ^^^2,n/ s y2{im-s) / s ^(s^2m) / s ^ ^ ^ (73) 

C = 70; [37 + 2 (a + 3/^ - 1)] + 2 [72 - 7 - 3/? (a + /?)] . (74) 
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• It turns out that C vanishes identicahy for aU solutions discussed in Cases 3 and 4. Therefore in these cases, 
integrating (I73p and setting the constant of integration equal to zero we find 

VW^i (75) 

-4to (w + 1) £171^2^^7-20) ^^^^ ^ ^ 2. 

Thus, the effective potential vanishes cither in the particular case where w = — 1 or in the Gaussian normal coordinate 
system (m = 0), i.e., for the solution ([5^ - (IMl) . 

• In Case 5, C 7^ and from ([75]) we get 

= <; ' (76) 

6/(2cj+ 1 - 2OI>A-2?;<,''^2"hn0, s = 2. 



• The parameter s is equal to zero only in Cases 1 and 2. In Case 1, the general- relativistic 5D solution (|46|) 
reduces on every hypersurface T,y to the spatially flat FRW models of general relativity, viz., 

with ^ 

P = Po ( T ) ' V = np, (78) 



(79) 



where po = 3/^/ (STrtg^l^yg) is the energy density measured at some time t — to and 

2-31 



31 



• In Case 2, to evaluate V we substitute (|49| into (IT9t . We find that (t)dV/d(f> vanishes identically by virtue of 
the equation relating /3 and w. The same is true for solution ([50]) . Therefore, without loss of generality we can set 
V = Q. Now, from ([72l) we find that (|49|) reduces on Ej, to a radiation dominated universe (p = 3p), regardless of w. 
Similarly, (|50l) reduces to empty space (p — p = 0) on Ey. 

4.1 4D cosmologies induced by the 5D solutions of Cases 3 and 4 

Cases 3 and 4, with s ^ 2, give more general BD cosmological model^ in AD. In fact, from ([7^ . using ([75)) we get 

8^ (BD) ^ m[(m-l)(s-2)-g2 + 1)] 
° (s-2)y2"^2i2 

ri(^^) = -nT0(^^\ (80) 



where 

s2 (w + 1) - Z (s - 2) 



(m-l)(s-2)-s2 (w + l)' 
Thus, the effective pressure and density satisfy the barotropic equation of state 



(81) 



p = np. (82) 



^If s = 2 and m ^ 0, the logarithmic potential leads to a violation of the weak energy condition p > 0. Therefore, to avoid such a 
violation we should set m = when s = 2, which in turns implies that T^^^^ = on every Sy. 
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We immediately notice that the 51? solutions (ISnt - lpO)) in the Gaussian normal frame (m — 0) reduce to scalar-vacuum 
cosmologies in 4£>. One can verify that for all solutions in Cases 3 and 4, n can be written as 



2 - 3Z - s 

n = , 

3/ 

and 



(83) 



m = 1 - 3? - s = -1 + 3Zn. (84) 
We can use these expressions to eliminate s and m in (|80p -(l8i p . viz., 

6/ (2 + 3n) - 6/2 [3n (n + 1) + 2] - 4 
[3/ (n + 1) - 2f ' 



Stt o(i3D) _ (/ + 2 - •jln) (1 - Zln) ^{bd) _ rr^oiBD) , . 

- (n + l)^2y2.^2 ' A - '^Jo ■ (85) 



As a consequence, we can formulate the reduced Brans-Dicke cosmological models in 4D, derived from the 5D 
solutions of Cases 3 and 4, solely in terms of I and n as 

_£ 

t 

^ X [2-3i(«+l)] 

3/(ri+l)/[3/(ri+l)-2] 

V = yo(:^j , (86) 



with 



- , to = ^J/o 



(; + 2-3;n)(l -3;n)^ 



16^po [?(l + 3n2)~2n] 
^° = 2 + Ul-3n) ' 



where po refers to the value of the energy density at some arbitrary fixed time to, and 



_ 6 (cj + 1) + 3n (2a; + 3) ± y/9n^ + 12 (1 + a;) (3n - 1) 
' ^^1^ 9(tJ + 2)n2 + 18(w + l)n + 3(3cj + 4) ' ^ ^ 

We recall that, in principle, I is an observable quantity because it is related to the acceleration parameter q = —ad/a^ 

by 

1 = ^. 

5 + 1 

The denominator of cj in vanishes when I = 2/3 (n + 1) and ut — >■ sign (3n — 1) x 00. In this limit s = 0, 
(j> = constant and from ((84)) it follows that m — {n ~ 1) / {n + 1), which means that the extra dimension contracts, 
or compactifies (m < 0), for any — 1 < n < 1. Consequently, formally for — >■ cx) we recover the usual spatially flat 
FRW cosmology of ordinary general relativity. 
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Matter-dominated universe: Since the present epoch of the universe is matter dominated let us consider the 
case where n = 0. In this case the weak energy condition is satisfied for any I > 0. On the other hand, the condition 

w > -3/2 restricts the range of / to be either I < 2 (^1 - y^2/tj « 0.37 or / > 2 (^1 + y^2/tj « 3.633. The former 
range requires q > (l + V6) /2 « 1.72, which is inapplicable to the present epoch. However, the latter range leads to 
q < (l — V6) /2 ~ —0.72, which is within the current observed measurements q — — 0.67± 0.25 [50]. We note that 
the parameter I is a real number only if (1 + < 0. 

Radiation-dominated universe: For n = 1/3 we find that Tq^^'^^ > in ([55]) requires I < 1, i.e., g > 0. Now 
from ([55)1 we get ^(^) = 1/2 and = 2 (1 + w) / (4a; + 5). The solution ~ 1/2 is identical to the radiation- 
dominated epoch of general relativity, regardless of the specific value of uj. The second solution violates the condition 
< Z < 1 for any uj < —1, as required by the dust model. 



Therefore, the effective BD cosmology in 4_D (|86 |) - (|88l) gives a decelerated radiation era, which is consistent with 
the big-bang nucleosynthesis scenario of general relativity (l = 1/2), as well as an accelerating matter dominated era 
compatible with present observations. 

To avoid misunderstanding, we should mention that the above discussion assumes that the parameters n and 
w are independent, which is true in the AD cosmologies generated by ([5^ and ([5^ : giving n and lu we obtain I 
from This in turn defines the rest of the solution ([55|) -(|57 1) . The 5D solutions ([5T|) and (157)) only have one 

independent parameter. Therefore they lead to a more restricted class of BD cosmologies in AD. In particular, one 
can show that (|5ip yields shrinking (instead of expanding) dust models. Regarding the solutions in the synchronous 
frame (|57p. they produce AD models of expanding dust, but do not yield accelerated expansion. 



• If we assume V = a.t the beginning, then we obtain a set of BD cosmologies in AD which precludes the 
existence of a dust- filled universe. In fact, in the case where V = 0, (m^0,a; = — 1) instead of ([50]) we have 

8n^Q{BD) ^ m {ni - 1) 

Thus, a dynamical universe (/ ^ 0) with p ^ requires p ^ 0. 

4.2 4D cosmologies induced by the 5D solutions of Case 5 

Finally, let us study the effective AD world generated by the 5D solutions discussed in Case 5. We emphasize that 
these are the only solutions demanding the extra dimension to be spacelike. 

Now the effective potential is given by (1761) . To avoid contradictions with the weak energy condition we disregard 
the logarithmic potential (see footnote 9). In this context, from ([7^ wc find the effective EMT as 

StT o(BD) _ 3/ [s {s + 31)) + (1 - s ~ 31)] liBD) _ rrOiBD) , . 

- ylt^^{s-2)A^ ' ^1 ' 

where n is the same as in Cases 3 and 4, namely 

2 - 3/ - s 



31 

As a consequence, from (|65l) we find that 7 can be expressed as 



^ ^ (3;n-2) [3;(n + l)-2] 

^ 3l{n + l)[3l{n + l){l+u)-Au:-3\ + 3l{l~\) + 2{l + 2uj)' ^ ' 
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Substituting this into (TOl) we get 



Sn oiBD) _ 2{2l~l- 2a;) + 3^ (n + 1) [4a; + 3 - ^ - 3^^ + 1) + 1)] 



After some calculations we find that the effective cosmology in 41? is formally identical to (|86|) except for the fact 
that now0 

. _ 87r(n+l)poil ,- _ . , 

to = Ayoto 



2 {21-1- 2uj) + 3/ (n + 1) [4cj + 3 - ; - 3; (n + 1) (a; + 1)] 



^ Snpo [31 (n + 1) - 2] {2 [1 + g; (1 - n)] + ^ [3a. (n^ - l) - 4] } 
° 2(2/-l-2a;) + 3/(n+l)[4w + 3-;-3;(n + l)(a; + l)] ■ ^ ' 

Let us consider the above sohition for large values of a;, in accordance with the solar system bound on a;(> 600). 
From the above expressions it follows that p > 0, > require / in the range 

2 2 1 

Consequently, as aj ^ oo we recover the usual spatially flat FRW cosmologies of general relativity. 

Let us now consider a matter-dominated universe (n = 0). Using that I = l/{q + l) we find that the positivity of 
the effective density requires 

9q -2q^ -1 



U! < 



i2q-lf 



For the sake of argument, let us consider an accelerated expansion with q = —2/3, which corresponds to Z = 3. This 
requires ui < —71/49 ^ —1.449. In general, q < requires uj < —1. As q — ^ 1/2, which is the general-relativistic 
value of q for dust, the theory becomes consistent with increasingly larger values of w. 

We note that the solution can give accelerating matter-dominated era as well as decelerating radiation-dominated 
era (n = 1/3) for the same value of ui in the range —3/2 < ui < —1. For example if we choose uj — —1.45, 
then a consistent model of decelerating radiation era is obtained for I < 0.88, and accelerated dust era for I > 3, 
approximately. 

• For completeness, let us consider the case where V = 0. Coming back to the original expression for (j>dV/d(j), 
we find that = requires 

^ {s + 3l-l)[s^ + s{3l-2) + 6lil-l)Y ^ ' 



Substituting into (j65p we obtain 

6(2/ - 1)/ 



s (s + 6/ - 2) " 

Thus, for V — (omitting intermediate calculations) we find 



(95) 



2-5-3/ 2(2/-l) 87r 2 - / + 37i/ [/ (3n + 1) - 3] 

p^np, n=^^, '^- (i_„)[2„3/(n+l)] ' T ° = v'AH^ (1 - n) " ^^^^ 

Consequently, the effective BD cosmology in AD can be written as 

^ \ 2(l-3n)/[4+3t<;(l-n^)] 

To) 

^ X 2[l+w(l-ri)]/[4+3Lj(l-n^)] 



«o ( - ) , (97) 



^'^Here ipo and Vb depend on three independent parameters, n, I and ui. If one fixes ui as in I I85II . then 1931 1 reduces to I I87II . 
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with 



47rpot^ [4 + 3w(l-n2)]' 



(3 + 2w) 4 - 6n + 3a; {n - ly 



Ay ok 



(98) 



In this parameterization, the n = models become identical to those presented in the original BD paper [38|. For 
n we recover the type A-I solutions discussed in 37 , although in a slightly different notation. From (^51) it is 
easy to see that any lu in the range —2 < uj < —3/2 (1 < / < 2) yields a matter-dominated universe with accelerated 
expansion. However, this range is inconsistent with a radiation dominated epoch, which demands lo — —3/2 (l — 1/2). 



5 Summary 

Brans-Dicke theory in AD can explain the observed accelerated expansion of the present matter-dominated universe, 
without invoking the presence of dark energy, if w is allowed to take some value in the range —2 < w < —3/2. 
However, uj in this range does not produce a consistent radiation-dominated epoch with decelerating expansion, 
as requires the big-bang nucleosynthesis scenario. One way out of this problem is to introduce a self interacting 
potential Another way is to consider a modified BD theory with a varying w |32| . Both approaches are not 

free of criticism. In one of them the potential is added by hand, while the other creates the necessity of finding the 
fundamental mechanism driving the variation of uj. 

A third way out of this problem is to resort to higher dimensions. In [40 it was shown that the BD field equations 
in 5-D can be reduced, on a hypersurface orthogonal to the extra dimension, to those of GR in AD coupled to two 
scalar fields. These fields may account for the late-time accelerated expansion provided several other conditions are 
met: the extra coordinate is cyclic; T^g ^ 0; x'^gh^ = 0; a'^$'" = constant, where m is a positive constant; the 
effective matter in AD has negligible pressure. 

In this work we have presented an approach which is free of the above criticism and conditions. Without making 
any assumption on the nature of the extra coordinate or the matter content in Si*, we have shown that the BD field 
equations in bD are equivalent to those of BD in AD derived from the action (I^Ul) with non- vanishing scalar potential 
V — V{<j)). The potential is not introduced by hand, instead the reduction procedure provides an expression, namely 
([T9l) . that determines the shape of V{4)) up to a constant of integration [ ([75l) - (|76|) ]. It also establishes the explicit 
formulae for the effective EMT in AD [ ([T4| -(fT6 |) ]. This extends and generalize some previous results recently obtained 
by the present author [53^. 

In the context of FRW cosmological models we have integrated the vacuum FE in bD under the sole assumption 
of separation of variables. We analyzed in detail the class of power- law solutions p3|) . We obtained a large family of 
solutions, namely (|67l) . that has three free parameters; various families with two free parameters [([50]). (15^ . ([5^ ]: 
one free parameter [ P^t . ([5T|) . ([57)) . ([5^ ]. and solutions that only exist for w = — 4/3 [ (p^ . Thus, the spectrum 

of BD power-law solutions in bD is significantly larger than the one of GR in 513, which can have at most one free 
parameter (Case 1). 

We discussed the effective AD world generated by our solutions. We found that the theory yields power-law and 
logarithmic potentials, except for the solutions derived in Case 2 and those in the Gaussian normal frame where we 
can set V — 0, without loss of generality. Certainly, one can assume V — aX the outset in (fT9|) . in which case the 
theory in AD would be identical to the original BD [38j. However, such assumption imposes a strong constraint on 
the parameters of the solutions leading to restricted cosmological models in AD which have no, or little, relevance to 
the problem of cosmic acceleration [ ([55]) . 

We found that all BD models in AD can formally be expressed as in l|86|) with 0o and Vb given by ([87]) or (|93p. The 
models are given in terms of the parameters n, I and u. When at least two of them are independent, which occurs 
for a large number of solutions [ (|50p . (15^ . ([55)1 . ([57|] we can obtain models that for the same a; can give the present 
accelerated expansion and a decelerated radiation-dominated epoch as required by primordial nucleosynthesis. 

We should mention that as in the conventional AD Brans-Dicke theory, in our models the (effective) EMT obeys 
the ordinary conservation law D^T^'-^^-' = (the same as in Einstein's theory), which in the cosmological realm 
yields the usual equation of motion 

p + 3-(p + p)=0. (99) 
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Although, we do not require the introduction ol a self interacting potential in AD or matter in 5D, the range of 
w obtained in our work is consistent with the one obtained by a number of authors in the literature in the context of 
several cosmological models and different versions of BD theory in AD and 5D [3T]-[3T]: In all cases cosmic accelerated 
expansion requires a; < — 1. 

There is also evidence in the literature that small values of |a;| are required by inflationary models [51J and structure 
formation |52j in scalar-tensor theory. Certainly, these theoretical results contradict solar system experiments which 
impose the constraint w > 600. However, this solar system limit does not necessarily imply that the evolution of the 
universe is, at all scales, close to general relativity. The fact is that GR is poorly tested on cosmic scales ^29, and 
no experiment has been done to test BD in cosmological scale yet. Consequently, one cannot discard scalar-tensor 
scenarios of the sort discussed here and in [31]- [41], [51]- [52] on the basis that in cosmological scales oj does not meet 
solar system requirements. 

Thus, the problem of (in)compatibility between astronomical and cosmological requirements remains open. Our 
recent work [S3j indicates that a variable (homogeneous) uj cannot solve the problem. Perhaps, local inhomogeneities 
might give rise to high values of uj consistent with astronomical observations |32) . 



Appendix A: Integrating the field equations 

In section 3 we have seen that the scalar- vacuum Brans-Dicke FE for the cosmological metric (|23p give five partial 
differential equations to be solved for four unknowns, n{t,y), a(t,y), ^(t,y) and (j)(t,y). The assumption that the 
metric is separable increases the number of unknowns to six. The aim of this appendix is to show, by means of 
explicit integration, the consistency of this assumption with the field equations. 



I: First, we assume that ci 7^ 0. Then from ((33|) we get 

where Qq is a constant. Substituting this and (j34l) into the wave equation (|28p we obtain 
ciFWW - W 



C3FW + {c2- ci)WF 



eClC2CsN§ [C2 + 3 (Cl - C3)] p-lry2 ^ q 

(3w + 4) 



(A-1) 



(A-2) 



When Cl 7^ the assumption of separability requires k = and the fulfillment of one of the following conditions (i) 
C — 0, (ii) 770 = 0, (iii) C2 — (cases la and lb bellow). 



la: Now we substitute the metric functions ([29]), including ([32]), ([34]) and (|A-ip . into the FE ([24])- ([26]) with fc = 
and either C = or 770 = 0. We find that ([M]) is a second-order differential equation for W, while (I^SI) and I^E^ 
contain the second derivatives of W and F. Next, we isolate W from (|A-2[) and substitute it into (|24l) to obtain the 
following first-order differential equation 

ciW^ , C2WF , ^ ^^^^^ 



with 



WF F^ F^ 



Cl = 3ciW + 6C1C3 — 2C3, 



3' 

C2 = 6C1 (C2 + C3 - Cl) - 4C2C3, 

63 = 2c2 (3ci - C2) , 

C4 = SecjCsN^ [m/ (3w + 4)) - 2cl] . (A-4) 

A similar substitution into (|26p yields an expression for F, whose explicit form we omit, a.s F = F{W, F,W , F). 
With these W and F we find that ([^5]) reduces to (jA-3p . As a result there are two independent equations, viz., (IA-2[) 
and (jA-3p . for the two unknowns W and F. 
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lb: For ci 7^ 0, C2 = (C 7^ 0,770 7^ 0), equation (jA-2p can be easily integrated. The general solution is = 
FoW~^W'''>^''^ , where Fq is a constant of integration. Using this expression, and (jA-ll) with C2 — 0, we find that 
both and (pS)) generate the same second-order differential equation for W, namely. 



1^2 

ci (ci - C3) M^PF + — 



CiO; + -c, — 



F2(3a. + 4) 



= 0. 



(A-5) 



On the other hand, ([^5]) becomes a third-order differential equation for W, which upon substitution of (jA-SP into it 
reduces to the identity = 0. 

• A particular solution to (jA-2p . for C2 = 0, is = constant, which leads to Q = constant. In turn, (|24l) and 
([26]) become identical to each other, both requiring C3 — ci, in agreement with (|A-5p . With this requirement (|25l) 
yields F oc t. 



II: We now consider the case where ci = and k ^ with C2 7^ and C2 = (cases Ila and lib bellow). 

Ila: First we assume C2 7^ 0, in which case from (1551) we get F — FqVF"^^/^^. Substituting this and ([M)) into the 
wave equation p8| we obtain^ ^1 

C2QWW - W UqW - ?,C2WQ) -f ^^ig^^g(^2- 3c3) 2(C.+C3)/C2 ^ 0^ (A-6) 

V / F^ (3w H- 4) 

When C2 7^ separability ([55)1 requires 770C = 0. Following the same steps as above, we obtain a first-order differential 
equation by substituting W from (|A-6p into Namely, 

Taking the time derivative of (jA-7p and replacing into it the expression for W obtained from (jA-6p . after some 
lengthly algebraic manipulations, we derive (j26p . which now is a second-order differential equation for Q. On the 
other hand, (l^5t reduces to (|A-7p after we substitute W and Q into it. Thus, (|A-6P - (|A-7p generate all the separable 
solutions with ci = 0, fc 7^ and C2 7^ 0. 



lib: We now let C2 = and assume C3 7^ 0. In this case one can generate solutions with A; 7^ and j/qC 7^ 0. From 
((33)) it follows that W = constant. Consequently, the wave equation (j28p is satisfied identically. From ((24)) and ((26)) 
we obtain the differential equations that govern F and Q. Namely, 

FQ , ^TVoV fc , ecl{Cs + C)N^ 



Q^ + FQ + li^J 2(3a.-K4)F2 " (^"^^ 

and 

Taking the time derivative of ([A-81) and using ((A-9p we obtain an expression for F, which we omit here, as F = 
Q{F,Q,F,W). Once again after some manipulations, one can verify that ((25)) reduces to ((A-Sp after F and Q are 
substituted into it. 

Thus, we have obtained the differential equations that generate all distinct cases of separable solutions. They 
demonstrate, in a constructive way, that the FE are consistent with the assumption of separability. In the case where 
ci = C2 = C3 = the metric functions become independent of the extra coordinate. For this case, the solutions to 
the FE as well as their cosmological applications have recently been discussed by the present author in [53] . 



'^^In general, when C2 7^ from II33I I we get F = FqW '^^^'^^Q 3ci/c2 g^^j ^^^^ wave equation II28I I yields 

C2QWW - W \c3QW + 3 (ci - C2) WQ] + ^^i'^^^o [^2 + 3(ci-C3)] ^(6,^+,^)/,^^2(,^+,3)/,^ 
L J {3lu + 4) 
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